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Abstract 

A hopping model for molecular motors is presented consisting of a state with 
asymmetric hopping rates with period 2 and a state with uniform hopping 
rates. State changes lead to a stationary unidirectional current of a particle. 
The current is explicitly calculated as a function of the rate of state changes, 
including also an external bias field. The Einstein relation between the linear 
mobility of the particle and its diffusion coefficient is investigated. The power 
input into the system is derived, as well as the power output resulting from 
the work performed against the bias field. The efficiency of this model is 
found to be rather small. 
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I. INTRODUCTION 



The modelling of molecular motors found great interest recently The aim is the 

understanding of the working of molecular motors, where motor proteins perform unidirec- 
tional motion on structures having potentials without inversion symmetry. Several processes 
were proposed that may lead to unidirectional motion of particles in such potentials. One 
process is a switching of the potential between different states, which may be caused by 
chemical reactions @,|10|. To understand the working of molecular motors in principle. 



simplified models are desired which are amenable to explicit calculations. In this note we 
present a very simple hopping model for unidirectional motion, which includes switching of 
the hopping potential between two states. Due to its simplicity, all quantities of interest can 
be calculated explicitly. 

Recently two other simplified hopping models for molecular motors were introduced. 
Schimansky-Geier et al. |TT[ considered chains of segments that can switch either individually 



or collectively between two states. The segments consist of three sites with a non-symmetric 
arrangement of hopping rates in one state and uniform hopping rates in the other state. 
Kolomeisky and Widom [|T2| studied a model with coherent switching between two periodic 
non-symmetric potentials. In each potential only downward transitions are allowed; con- 
sequently the motion of a particle comes to rest in each state in the absence of switching 
events. Our model is simpler than the model of Schimansky-Geier in that the segments have 
only length 2. We restrict the derivations to coherent switching between two states. The 
model is also different from the one of Kolomeisky and Widom in that it includes diffusion 
of the particle in the two states; unidirectional motion is then brought about by the state 
changes. 

We restrict in this note the derivations to the stationary state with a unidirectional 
particle current; also we do not intend to model real molecular motors with our simplified 
model. In the following section the model is formulated in terms of its master equations. In 
section III the stationary current is calculated. The case with an additional uniform bias is 
treated in Sec. IV. The question of the validity of the Einstein relation between mobility and 
diffusion coefficient is addressed in Sec. V. Section VI presents a derivation of the power 
input and the power output when the systems performs work against the uniform bias field. 
Section VII contains concluding remarks. 



II. DEFINITION OF THE MODEL 

The model consists of a linear chain with sites to which hopping rates to the neighbor 
sites are associated. The chain exists in two different states with different hopping rates. The 
model is pictorially represented in Fig.l where the hopping rates are visualized by potentials. 
In the state (1) with an asymmetric potential, the sites with rates 71, 5i are designated by 
1, and the ones with 72,^2 by 2. The transition to the other state (2) occurs in form of a 
Poisson process with rate Tf, and a particle on a site of type 1 or 2 is transferred to a site 
designated by 3 or 4, respectively. In the second state with uniform potential all sites have 
transition rates 7 to the neighbor sites. The reverse transition to the non-symmetric state 
occurs with the same rate Tf. Periodic boundary conditions will be imposed on the system. 
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To derive the current, it is sufficient to take a period 2. Displacements of the particle, 
however, have to be counted in the extended potentials. 

Let Pi{t) be the probability of finding the particle in any of the sites of type 1 in the 
chain in state (1). The probabilities ^2(^)5 Psit) and P^it) are defined analogously. The 
probabilities fulffil the master equations 

(72 + P2 - (71 + Si) Pi + Tf{Ps- Pi) 
(71 + Si)Pi - (72 + ^2)^2 + TfiP, - P2) 
27P4-27P3 + r/(Pi-P3) 

27P3-27P4 + r;(P2-P4) (1) 

The initial conditions are of no interest for the stationary state. 

The sum of the four equation shows that the time derivative of the sum of the four 
probabilities vanishes. Hence Pi + P2 + P3 + P4 = const = 1. The sum of the ffist and the 
second master equation, and the sum of the third and fourth equation lead to the system of 
equations 

r/(P3 + P4)-rKPi + P2) 

r/(Pi + P2)-r^(P3 + P4). (2) 

This set of equations shows that the state changes occur with rate Tf. The stationary 
solution of these equations is P1+P2 = P3+P4 = | where we have also used the normalization 
of the sum of all probabilities. 



dPi 

~dt 
dP2 
dt 

dt 
dPi 

~~dt 



djPi + P2) 
dt 

djPs + Pa) 
dt 



III. STATIONARY CURRENT 



The current from sites of type 1 to 2, or from 3 to 4, respectively, in the positive direction 
is given by 

J21 = 71P1 - 52P2 + 7(^3 - Pa) . (3) 

The current from 2 to 1, or from 4 to 3, respectively, in the positive direction is given by 

J12 = I2P2 - SiPi + -f{PA - P3) . (4) 

In the stationary situation, the two currents are equal, J21 = J12 = J- Equating Eqs. (j^) and 
(1) and using the stationary values for the sums Pi + P2 and P3 + P4, we obtain an equation 
relating the differences P2 — Pi and P4 — P3, 

(71 + 72 + '^1 + ^2)^(^2 - Pi) + 27(^4 - P3) = \ili -12 + Si- 52) . (5) 
We find a second equation by subtracting the fourth from the third equation of the set (P, 
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Tf{P2 - Pi) - (47 + r;)(P4 - Ps) = . (6) 

The set of the two equations and for the two differences is readily solved. We proceed 
directly to the stationary current, which follows from the sum of Eqs. (^) and (^) by inserting 
the differences, and the stationary values for the sums. The final result for the current is 

J ^ 1 7r/(7i + 72-^1 -S2) 

2 (71 +72 + ^1 + 52)(47 + r/) +47!/ ■ 

In the derivation of the current, there appears a term in the numerator proportional to 
7i72 — ^1^2, see also Eq.@ below. This term disappears due to the condition of detailed 
balance, which is valid within each of the two states of the chain. Consequently, the current 
is proportional to for small Tf, as it should be. We point out that the condition of 
detailed balance is fulfilled within each of the two states. The introduction of the state 
change with the rate Tf leads to a violation of the condition of detailed balance when the 
full system is considered. This violation is necessary for the working of the motor, cf. the 



discussion in Ref. |jTO|. 



The result shows that a unidirectional current exists for finite values of the flip rate Tj. 
Figure 2 presents the result for a particular choice of the hopping rates on the chain, as a 
function of the flip rate Tf. The figure also contains the results of numerical simulations of 
the model. The current saturates for Fj >> 7 to the value 

J 1 7i + 72 - ^1 - S2 , . 

2 ^ 7i + 72 + 5i + 5, + 47 • 

From the continuum models with asymmetric potentials one expects that the unidirectional 
current is maximal for an optimal flip rate and vanishes in the limit of infinite flip rate. 
The saturation of the current at infinite flip rate is an artifact of our hopping model with 
only two sites. The model would correspond, in a continuum version, to a model containing 
sections with infinite slope in the potential. Models with at least three sites are required 
to remove the artifact of saturation. For instance, the model of Schimansky-Geier et al. 



11| exhibits a maximum of the current as a function of the flip rate. The vanishing of the 
current in the limit of infinite flip rate, however, requires the presence of infinitely many 
sites in a segment, i.e., the continuum limit. 



IV. UNIFORM BIAS 

It is also instructive to consider the influence of a uniform bias on the unidirectional 
motion of the particle. Then the system has to work against the bias field and the question 
is which bias is sufficient to stop the motion completely. A bias is introduced by multiplying 
the transition rates in state (1) in positive direction by a factor 6 > 1, 71,72 — * &7i,&72, 
and the transition rates in negative direction by b~^, 61,62 —>■ b^^6i,b~^62- In the state (2), 
the transition rates 7 are multiplied by b for the positive direction and b~^ for the negative 
direction. The bias factor is related to an external force F hy b = exp[Fa/2kBT) where a 
is the lattice constant which is usually set unity. 

The calculation of the stationary current is easily extended to the case with bias, using 
the modified transition rates as given above. The result is 
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+ (6 - 6-i)7[&(7i + 72) + h-\5r + ^2)]} 

+ (6 + ri)7r^[6(7i + 72) - ri(5i + 52) + 2(6 - 6-^)7]) (9) 



with 

A(6) = [6(71 + 72) + h-\6i + 52)] [2(6 + 6-1)7 + T^] + 2(6 + h-^)^T j . (10) 

The resuh (|^) is readily recovered for 6=1. 

The result for the stationary current as a function of the bias 6 is presented in Fig. 3, 
for one particular value of the flip rate Tf. One notices that the current behaves almost 
linearly as a function of 6, for our choice of the parameters. The bias value where the current 
vanishes can be read from the figure. We note that this value is rather small. No analytic 
expression can be given because the condition J(6) = corresponds to an equation of G*'* 
order in 6. The slope of the curve is given by 

% = A(6^T) t^^i^' + '^i'^2)(47 + T/) + 7(27 + r/)(7i + 72 + 5i + ^2) + 27'r^] 

~ 2A(6 = 1)2 [^^^(^^ + ^f^^^' + 72 - ^1 - ^2)'] (11) 

Using the slope and the value of the current at 6 = 1, an approximate prediction for the 
value of 6 where the current vanishes can be obtained. 



V. MOBILITY AND DIFFUSIVITY 

The slope of J(6) at 6 = 1 gives a prediction for the linear mobility of a particle in the 
stationary situation, i.e., the increase/decrease of its velocity due to the action of a driving 
force. The question arises whether the linear mobility B can be related to a diffusion 
coefficient D by the Einstein relation D = ksTB. 

To examine this question, we calculated the diffusion coefficient of a particle in the 
system from its asymptotic mean-square displacement where the squared displacement is 
subtracted, 

<x^{t)> - <x{t)>'^ — ^2Dt 7t->oo. (12) 

To derive the mean-square displacement, site probabilities have to be considered that are 
defined on a linear chain with 2N sites, in the limit N —>■ 00. The master equations (|l|) are 
easily extended to include site indices. The derivation of the mean-square displacement is a 



standard procedure, see for instance []T3| or Because it is somewhat lengthy, it will not 
be reproduced here. 

The result for the asymptotic diffusion coefficient is 

D = ^(^^[(7i72 + 5i52)(47 + Tf)+ 7(27 + r^)(7i + 72 + + ^2) + 2^''T j] (13) 
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The diffusion coefficient is identical to the first term of the hnear mobihty as given in 
Eg ■ (pTl) ■ Since the second term is missing, the Einstein relation between mobility and 
diffusion coefficient does not hold here (note that we are using units where ksT = 1). A 
violation of the Einstein relation is not surprising for the driven non-equilibrium, although 
stationary, system that is considered here. We point out that the numerical difference 
between the mobility and the diffusion coefficient is rather small for the parameters chosen 
here, it is barely visible on the scales of Figures 3 or 4. 

Two limiting cases of the diffusion coefficient are of interest. First, for small flip rates D 
can be brought into the form 

^^^f ^ 0) = ? + ^-1 ^ ^-1 1 ,-1 ,-1 • (14) 

This form of the diffusion coefficient follows from a two-state model where D is deduced 
from a weighted superposition of the diffusion coefficients in state (2) (first term on the rhs 
of ([T^)) and of the coefficient in state (1) (second term). Since the flip rate is the same in 
both directions, the weights are 1/2, and the result that is obtained agrees with (|^) and 
is independent of the flip rate. Our derivations show that a two-state description is only 
applicable at small flip rates. For large flip rates, the diffusion coefficient saturates at 

D{Tf^ 00) = 'J z — . (15) 

71 + 72 + 51 + ^2 + 47 ^ ^ 

We were not able to find a physical interpretation of this limiting behavior. 

We also performed Monte Carlo simulations of the mean-square displacement of particles 
about the stationary state. The number of particles and system size used in the simulation 
are respectively 10^ and 4 with periodic boundary condition. Figure 4 demonstrates that 
< x'^it) > — < x(t) >^ is indeed proportional to t , and the diffusion coefficient D can be 
deduced from the simulation results. The numerical results agree with the theoretical ones. 



VI. POWER INPUT AND OUTPUT 

In this section we calculate the power input when the system is switched between two 
states, with an additional bias acting on the particle, and the power output, i.e., the work 
performed by the particle per unit time when it moves against the uniform bias. 

Power is required to switch the system from the state with a uniform potential to the 
state with the non-symmetric potential. In our model, if a particle is on a site of type 4 
in the uniform potential at the time where the switching occurs, its energy will be raised 
when it finds itself on site type 2 after the switching event. The energy follows from detailed 
balance in state (1), 

/3e = ln^ = ln^. (16) 

7i "^1 

To obtain the power input, the average occupation of sites type 4 at the switching 
event is required. To derive this quantity, we consider the dynamics of the site occupation 
under the condition that the particle is remaining in state (1), or is remaining in state 
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(2). The respective conditional probabilities, which will be called Qi{t),Q2{t) for state (1) 
and Qsit), Qiit) for state (2), obey analogous master equations as Pi(t), . . . , P4{t), cf. (P, 
except that the terms proportional to are absent. The two uncoupled systems of master 
equations are readily solved and the solutions are 

Qi(t) = A, + Bi exp[-(7i + 72 + 5i + ^2)^] 
Q2{t) = A2 + B2 exp[-(7i + 72 + 5i + ^2)^] 
Q^{t) = As + Bsexp{-A^t) 

Qiit) = Ai + Bi exp(-47t) . (17) 



The constants are given as 



72 + §2 



7i + 72 + (^l + ^2 
(7i + '^i)Qi(0)-(72 + ^2)g2(0) 
7i + 72 + ^1 + S2 
7i + ^1 
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Ay 

7i + 72 + c)i + 62 
B2 = —Bi 

As = A4 = ^ 

B; = -B, = ^[Qs{0)-Q,m. (18) 

The strategy to obtain the average occupation probabilities of the four types of sites at the 
average switching times is as follows: We use as initial conditions for Qi{0) in the solutions 
of the uncoupled systems of master equations the complementary occupation probabilities 
Qjitf) in the other states at the mean switching time t/, e.g., Qi(0) = Q^{tf). The mean 
switching time is given by = l/Fj. The coefficients for the solutions of the uncoupled 
master equations are then 

A, 
Bi 



72 + ^2 



7i + 72 + ^1 + ^2 

{li + 5i)Q-i{tf)-{i2 + 52)Q4{tf) 

7i + 72 + (^1 + ^2 

7i + ^1 



2 



A^- r- r- 

7i + 72 + ^1 + o- 
B2 = —Bi 

1 

A3 = Ai = - 

B, = -B, = hQ,itf)-Q2itj)]. (19) 



We now consider the solutions Eq.(|T^ at the time tj using the coefficients given in Eq. (p^) . 
They constitute sets of hnear equations for the required coefficients Qj{tf). The solutions 
of these equations are 
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where 
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+ 47, 


(21) 



We have verified by numerical simulations, where the system was switched at random times 
according to a Poisson process with mean fiip time tf, that the expressions given above 
describe the average occupations at the switching times. The above derivations were formally 
made for the case of no bias. A uniform bias is immediately included in the derivations if the 
transition rates 71, 72 are replaced by 671, 672, and the transition rates 61, 62 by b~^6i, b~^62, 
respectively. The final formulae will be considered for the case including bias. 

The power input is now obtained by multiplying the occupation of type 4 sites at the 
fiipping time with the fiipping rate and the energy supplied in one fiip. This energy has 
been given in Eq.(0) in units of fc^T. Hence the power input is, in units of inverse time 

Pin = Q4itf)Tfe/2kBT. (22) 

The power output is obtained from the velocity induced by the switching process by 
multiplying it with the force. The force is obtained from the bias parameter as 

F/keT = 2 ln(6) (23) 

where the expression is dimensionless due to the division with /c^T and the use of lattice 
constant a = 1. Since the velocity is related to the current by a factor 2 in view of the 
normalization of 1 particle in the reduced system with periodic boundary conditions, we 
have 

Pout = ±2J{F) F/ksT . (24) 

If the current for F ^ is negative for the choice of parameters used (as it is the case here) 
we compensate the sign by multiplying with —1. The current is obtained as a function of F 
from Eq.(H) by using Eq.(^). 

The power input, the power output and the efficiency of the motor, which is usually 
defined as the ratio of the power output to the input, are plotted in Fig. 5, 6 and 7 respectively 
for three chosen fiip rates as a function of the force acting on the system. The increment 
of the power input is so small that it cannot be seen from the figure. For example, for 
Tf = 0.2 the power input increases from 0.199 993 36 to 0.199 994 01 in the range of the 
force considered, from to 0.1. The efficiency of the model that is studied here is rather 
small. This is related to the fact that the unidirectional current vanishes already for small 
bias fields opposing the motion of the particles, cf. Figs. 3 and 6. In order to devise models 
that have a higher efficiency, more sites within a repeat unit are required, as well as proper 
adjustment of the parameters. 
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VII. CONCLUDING REMARKS 



In this paper we have presented a very simple model for a molecular motor where a 
particle moves unidirectionally in a hopping potential which is flipping between a state 
having a potential without inversion symmetry and a state with uniform potential. Such 
models are intended to mimick motion of proteins in polar filaments [|10| . The appeal of the 
model, which comprises only two different types of sites in the non symmetric state, is that 
all quantities of interest can be calculated explicitly. We demonstrated this by calculating 
the stationary current with and without bias which results from the switching of the system 
between the two states, and by deriving the power input and output of the system. The 
efficiency of the system as a motor is rather small for the parameters that were used. It 
was not the intention in any case to model real molecular motors. We also compared the 
linear mobility under the influence of a small bias field with the diffusion coefficient that 
follows from the mean-square displacement. We observed a violation of the Einstein relation 
between mobility and diffusion coefficient. We point out that we have a stationary non- 
equilibrium system where the question of the validity of the Einstein relation is nontrivial. 
The investigation of the fluctuations in these non-equilibrium systems appears to be of 
further interest. 
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Figure Captions 



Fig. 1: Pictorial representation of the model for a molecular motor, with the transition rates 
indicated. 

Fig. 2: Particle current without bias as a function of the flip rate P/ between the two states. 
Line: analytical result, symbols: simulation results. The value of the transition rates used 
in the simulations are: 71 = 72 = 7 exp(— 2), 5i = 7 exp(— 4) , 7 = ^2 = 0.5. Transition 
from state (1) to state (2) and vice versa occurs with the rate P/. 

Fig. 3: Particle current as a function of the bias factor b for a fixed fiip rate, P/ = 0.5. Line: 
analytical result, symbols: simulation results. 

Fig. 4: Mean square displacement, where the squared mean displacement is subtracted, 
shown for different flip rates as a function of time. Symbols:* simulation results for P/ = 5.0; 
X simulation results for P/ = 0.5; + simulation results for Tf = 0.1; short dashes: theory 
for Tf = 5.0; long dashes: theory for Tf = 0.5; full line: theory for Tj = 0.1. 

Fig. 5: Results for the power input for three selected values of P/. The corresponding values 
of the flip rates are indicated in the figure. Force is measured in units of ksT 

Fig. 6: Power output for three different values of P/ in units of ks T. The values of P/ are 
indicated in the figure. 

Fig. 7: Efficiency of the system vs force for three selected values of P/, which are indicated 
on the curves. 
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